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ABSTRACT 


A  new  and  simple  device,  similar  to  an  abacus,  was  invented  by  a  Chinese 
mathematician;  it  is  used  for  multiplication,  division,  and  extracting  square  roots. 
The  underlying  idea  is  abstracted  and  generalized  for  possible  application  to  con¬ 
temporary  computers  in  order  to  reduce  computing  time. 
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SECTION  I 


INTRODUCTION 

A  very  simple  device  was  recently  invented  in  China  for  effecting  numerical 
multiplication,  division,  and  extracting  square  roots,  all  without  the  use  of  the 
multiplication  table. 

According  to  a  popular  journal  [Reference  -cm]  from  which  this  information 
was  obtained,  the  inventor,  the  mathematician  YU  Chan-Hsan,  was  originally  a 
carpenter  and  later  a  clerk  in  a  stationary  store.  A  movie  of  his  biography  was 
produced  because  his  invention  is  comparable  only  to  that  of  the  abacus.  His  bar- 
multiplier  is  to  the  process  of  multiplication  what  the  abacus  is  to  addition,  and 
just  as  the  abacus  is  universally  used  in  China,  YU’s  bar-multiplier  in  the  combined 
form  with  the  abacus  is  also  becoming  part  of  the  standard  tools  of  his  country. 
Because  it  is  so  simple  and  easy  to  learn,  it  is  now  taught  at  almost  all  grammar 
schools  of  Mainland  China. 

TTNapierTs  Bones, n  invented  by  John  Napier  (1550-1617),  originator  of  numerical 
calculation  by  the  logarithm,  are  basically  identical  with  the  method  of  Figure  1 
[References  -EB:cm,  N-WM1,  46l].  The  new  Chinese  method  differs  from  it  in 


A  =  1 

4 

2 

3 

X 

B  = 

8 

3 

2 

2 

8 

4 

6 

A.  B. 
i  J 

+  8 

3 

32 

12 

16 

6 

24 

9 

8 

35 

30 

38 

13 

6 

3 

5 

3 

8 

6 

8 

3 

0 

1 

3 

V 

II 

H 

1 

8 

3 

9 

3 

6 

Figure  I  An  Example  of  Multiplication 
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the  use  of  the  bars  instead  of  the  numerals  and  in  combining  the  bar-multiplier  with 
the  abacus  as  the  " bar-bead  calculator."  It  converts  the  arithmetic  operations  into 
discrete,  purely  geometric  operations  on  points  and  lines  such  as  embodying  lines, 
by  determining  intersections  in  the  bar  part  and  linear  displacement  of  beads  in  the 
abacus  part. 

This  document  explains  YU!s  bar-calculator  as  a  prototype  of  the  new  idea,  and 
then  generalizes  it  to  the  multiplication  of  many  numbers  so  the  numerical  computa¬ 
tion  of  the  N-th  root  becomes  simple  and  fast. 

Application  of  the  basic  principle  of  the  prototype  gives  a  possibility  of  reducing 
computing  time  with  some  increase  in  computing  elements.  For  the  generalized  ease, 
further  reduction  is  possible  by  a  further  increase  of  the  computing  elements.  The 
key  lies  in  the  separation  of  the  pure  multiplying  operation  of  individual  digits  from 
the  carrying-over  process,  whieh  can  be  performed  in  the  subsequent  adding  step. 

When  an  equation  is  valid  only  for  the  multiplication  of  N  =  2  or  3  digits,  namely 
AB  or  ABC,  it  is  shown  by  adding  the  denominator  N  to  the  number  of  the  equation 
whieh  is  valid  for  the  general  value  N.  Variables  are  indicated  by  lowrer-case 
letters;  constants  are  indicated  by  eapital  letters.  All  lower-ease  letter  indices  arc 
variable  and  all  upper-case  letter  indices  are  fixed. 

From  the  standpoint  of  modern  abstract  algebra,  if  the  ordinary  earried-over 
produet  P  is  substituted  by  the  "quasi-deeimal,"  "quasi-binary"  or  any  "quasi" 
digital  expression  Q  of  the  product  defined  below  by  withholding  the  carryover  proc¬ 
ess,  a  lueid  isomorphism  between  the  digital  numbers  of  any  single  base  and  the 
functions  of  a  single  variable  is  revealed. 

As  a  whole,  this  paper  means  the  realization  of  direct  isomorphic  mapping  of 
the  arithmetic  operations  on  numbers  into  discrete  affine  geometric  operations  sueh 
as  meet  or  join  on  points,  lines,  planes,  or  hyperplanes  by  means  of  the  logical 
electrie  circuits. 

Though  some  literature  [C-GN,  F-SZ,  H-MGN,  L-DG,  M-GZ]  on  the 
"geometry  of  numbers"  already  exists,  none  includes  the  geometrical  method 
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of  multiplying  digital  numbers  by  discrete  models,  whieh  also  seems  to  deserve 
the  term  "geometry  of  numbers,"  Closely  related  problems  are  diseussed  in 
computer  design,  but  without  geometrie  interpretation  [it-ADC] . 

The  generalization  is  easily  obtained  by  considering  the  diagonal  hyperplanes 
expressed  by  the  eovectors  in  the  N-dimensional  affine  space.  Though  the  following 
explanation  is  given  mostly  on  the  deeimal  ease,  the  result  is  valid  for  any  base, 
and  would  be  particularly  advantageous  in  the  binary  system. 
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SECTION  II 


PRINCIPLE  OF  TIIE  BAR-MULTIPLIEIl 


For  the  multiplication  of  two  numbers: 


n  =  A,  B 
P  =  AB 


1/2 

2/2 


where 


S  s 

P  =  Po10  +  .  .  .  +  P  10  +  .  .  .  +  P  ,  P0  /  0,  P  <  10  , 
S  s  o  S  s 


s  =  0  ,  1  ,  ...  ,  S ; 


A  - 


=  AjlO  +  .  .  .  +  A.101  +  . . .  +  Aq  ,  Aj  /  0  ,  A.  <  10  , 


5.  1 


i  =  0  ,  1  ,  I  ; 


6.  1 


B  =  B  10J  +  .  .  .  +  B.10  +  .  .  .  +  B  ,  B  /  0  ,  B.  <  10  , 
J  j  o  J  j 


5.  2 


j  =  0  ,  1  ,  . .  .  ,  J  ; 


6.  2 


if  the  products  A.  B^  for  all  pairs  of  digits  A^  and  B^  are  written  down,  deferring 

the  carry-over  from  A.  B.  to  A.  ,B.;  e.  g.  ,  1  of  A„B,  =  4.  3  =  12  of  Figures  1  and 

i  J  i+l  J  21  B 

2toAB  =1.3=3  into  the  final  adding  step,  the  multiplying  process  will  change 
3  1 

to  the  form  shown  in  Figures  1  and  2. 


In  the  binary  system  there  is  no  carry-over  in  the  partial  product  AB^ .  In 
Figure  1,  the  two-digit  number  of  the  sum  of  each  column  are  vertically  staggered 
with  underlines  for  the  convenience  of  the  succeeding  carry-over  step.  Thus,  the 
products  of  the  type  A^B^  are  perfectly  separated  from  each  other,  and  are  mutually 
independent  in  logic.  By  this  separation,  the  geometrical  location  of  the  factors 
A  and  B  in  Figure  2  becomes  very  simple:  all  B  appear  in  the  first  row  of  the 

i  j  ° 
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Figure  2  The  General  Case 


uncarried  partial  product  and  all  appear  in  the  first  right  skew  line.  Therefore, 

as  shown  in  Figure  3,  if  two  lines  are  drawn  horizontally  through  the  positions  of  all 

and  three  skew  lines  are  drawn  through  the  positions  of  Aq  ,  then  the  number  of 

the  cross-points  of  these  lines  at  the  intersection  corresponding  to  A^B^  is  exactly 

A  B  =  2.  3  =  6.  Thus,  if  all  numerals  A.  and  B.  are  substituted  by  lines  as  shown 
o  o  1  ] 

in  Figure  3 ,  the  number  of  the  cross-points  at  all  of 


M  =  (I  +  1)  (J  +  1)  =  4.3  =  12 

Ld 


intersections  gives  the  values  of  the  products  A.B^  • 


If  both  are  D  digits: 


I+1=J+1=D, 


7/2 


8/2 


then  , 


M2 


9/2 
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Figure  3  Bar-Multiplier 

Because  A  and  B  are  commutative,  by  inclining  the  lines  of  B^  as  shown  in 
Figure  4,  A  and  B  are  brought  into  symmetry.  If  the  bar- multiplier  is  combined 
with  a  Japanese  abacus  (which  has  only  five  beads  for  each  digit  — one  upper  bead 
expressing  five  and  four  lower  beads  expressing  unity),  the  number  of  bars  for 
each  digit  can  be  similarly  reduced  to  five. 

Further,  instead  of  removable  bars,  flat  bars  can  be  installed  which  can  be 
turned  a  certain  angle  by  small  staggered  levers  so  that  each  digit  can  be  expressed 
by  an  angular  flipping  of  these  bars.  Now  based  on  Figure  4,  the  principle  of  the 
"bar-multiplier"  is  stated  as  follows: 

If  two  numbers  A  and  B  ,  e.  g.  ,  1423  and  832,  are  to  be  multiplied,  place  or 
draw  one,  four,  two  and  three  parallel  inclined  bars  from  the  left  lower  position 
to  the  upper  right,  and  eight,  three,  and  two  inclined  bars  from  the  left  upper  posi¬ 
tion  to  the  right  lower  as  shown  in  Figure  4  then  add  the  number  of  the  cross-points 
along  the  vertical  columns  with  carry  over,  to  obtain  the  product  P,  e.  g.  ,  1183936. 
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Figure  4  Bar-Multiplier  in  Symmetric  Form 


If  p  -  H83936  and  A  =  1423  are  given,  B  can  be  determined  by  subtracting 

the  cross-points  generated  by  the  maximum  digits  B..  Thus  division  is  possible 

as  exemplified  for  68373^-273  in  Figure  5.  Similarly,  square  and  square  root 

2 

can  be  performed  by  this  bar-multiplier  as  exemplified  for  213  in  Figure  6  and 
1/2 

for  45369  in  Figure  7. 

2P  PROOF 

In  the  products  A.B.101  +\  at  first  the  sum  of  all  the  terms  of  power  t 

(see  Figures  1  and  2)  is  defined  by 

Q  =2.  .  R..  ,  10/2 

t  i,  j  ij 
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Figure  5  Dividing  68373  by  213 


Figure  6  Squaring  213 


R..  =  A.B.  11/2 

1J  i  J 


over  all  values  of  i  and  j  of 

6.1,  6.  2:  i  =  0,  1,  ...  ,  I;  j  =  0,  1 . J  ; 

satisfying 


i  +  j  =  t  , 


12/2 


t  =  0,  1,  ...  ,  T  ; 


13 


T  =  I  +  J  . 


14/2 


Q  is  defined  as  the  sum  of  all  Q^IO  : 


Q  =  Q^IO^  +  . . .  +  Qt10  +  . . .  +  Qq 


15 


Though  Q  is  equal  to  the  product  P: 


Q  =  P  , 


Q  may  generally  be  not  only  larger  than  ten: 


Qt ; 10 


16 


17 


9 


Figure  7  Square  Root  of  45369 


as  seen  in  the  example  of  Figure  1,  but  it  may  consist  of  three  or  more  digits. 
Hence,  Q  will  be  called  the  ”  quasi-decimal”  expression  of  the  produet.  Now,  on 
the  plane  of  the  paper  of  Figure  4,  if  one  introduces  an  affine  frame  (see  the  end 
of  this  paragraph)  having  the  bases 

=  0A1,  "a2  =  0A2  a- 1/2 
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where  the  origin  0  is  at  the  point  R^,  the  first  base  point  A^  is  at  the  point  R 
the  second  base  point  A  is  at  the  point  R  (as  shown  in  Figure  8),  then  a  linear 

&  U  X 

equation 

1  2 

w  x  +w  x  =  1  a- 2/2 

on  the  components  of  a  vector  x  where 

—  1—2—  „ 

x  =  x  a  +  x  a2  a-3/ 2 

expresses  a  straight  line  intersecting  with  1-axis  at  l/w  and  2-axis  at  1/w  and 

X  ^ 

1  2  . 
w  x  +  w2x  =  t;t  =  0,  1 . T  a-4/  2 

express  straight  line  in  parallel  to  that  of  equation  a- 2/2  with  the  intercepts  t/w 
and  t/w'2  on  1-axis  and  2-axis,  w  and  w  are  covariant  referring  to  x*1  (n  =  1,  2) 
and  called  the  components  of  a  covector  w  : 

12  „ 

w  =  w  a  +w  a  a-5/2 

►  X  *  z  ^ 


where 


a-6.  1/2 


is  the  base  covector  expressed  by  a  pair  of  an  initial  straight  line  0  passing  through 

2 

the  origin  and  a  terminal  straight  line  A  passing  through  the  base  point  A  ,  both 
in  parallel  to  2-axis. 


Thus, 


a-G.  2/2 


is  the  base  covector  expressed  by  a  pair  of  an  initial  straight  line  0  passing  through 

2 

the  origin  and  a  terminal  straight  line  A  passing  through  the  base  point  A  ,  both 

12  2 

in  parallel  to  1-axis.  These  base  covectors  a  and  a  are  contravariant  as  ex¬ 
pressed  by  their  superscript.  If  the  components  of  the  covector  w  are  both  unity, 


11 


Figure  8  Affine  Coordinates  for  Figure  4 


12 


a- 7/ 2 


(w1,  w2)  =  (1,  1)  , 


the  corresponding  equation 


a- 8/ 2 


expresses  the  straight  line  passing  the  base  points  and  A0  ;  namely,  a  part, 
A  A  ,  of  it  forms  the  diagonal  of  the  base  parallelotope.  Similarly, 

X  z 


a- 9/ 2 


expresses  the  straight  line  passing  the  points  (t,  0)  and  (0,  t)  in  parallel  to  that  of 
equation  a-8/2. 


By  comparison  of  equations  12/2,  13,  and  a-9/2,  we  see  that  is  the  sum  of 


the  products  R. .  lying  on  the  diagonal  lines  of  equation  a-9/ 2  as  seen  in  Figures  4 
and  8. 

Thus  the  validity  of  the  above  principle  was  proved  by  introducing  the  affine 
frame.  If  this  affine  frame  is  used,  the  above  process  can  be  stated  as  follows  as 
the  case  2  of  the  general  case  N  described  below: 


1.  Construct  a  two-dimensional  affine  frame; 


2.  Enter  A^  to  all  points 


2 


=  j  ,  j  =  0,  . . .  ,  J 


18.  1/2 


x 


on  the  line 


x1  =  i  ,  i  =  0,  .  .  .  ,  I  ; 


19.  1/2 


3.  Enter  B.  to  all  points: 


1 

x  =i,i=0,...,I 


18.  2/2 


13 


on  the  line 


2 

x 


=  j  , 


j  =  0,  .  .  .  ,  J  ; 


19.  2/2 


4.  By  adding  all  products  R  on  the  diagonal  lines 
,  1  2 

a-9/2:  x  +x  =  t  ,  13:  t  =  0,  . .  .  ,  T  ;  14/2:  T  =  I  +  J  20/2 

and  hence  Q  are  obtained; 

5.  By  carrying,  P  is  obtained  from  Q. 

A  note  on  the  "affine"  space  [see  W-STM,  S-RC  1,  Wh-GI  349-53,  G-LA  250, 

A-GA  66]  :  If  different  lengths  are  used  for  units  of  different  axes  of  an  oblique 
coordinate  system  in  N-dimensional  Euclidean  space,  namely  a  "heterometric" 
oblique  frame  is  used,  the  frame  can  be  called  the  "affine"  frame  without  loss  of 
generality.  The  unit  of  length  of  any  direction  other  than  those  of  the  axes  is  not 
defined  and  is  perfectly  arbitrary.  The  inner  (or  dot)  product  is  defined  only  for  a 
vector  v  =  0V  and  a  covector  w  =  0W  where  0W  means  a  pair  of  initial  and  terminal 
hyperplanes  in  parallel,  with  intercepts  l/wn  on  n-axis  (n  =  1,  2,  .  .  .  ,  N)  (straight 
lines  for  N  =  2,  planes  for  N  =  3).  When  the  origin  is  fixed,  the  initial  hyperplane 
passing  through  the  origin  is  often  omitted.  The  value  of  the  inner  product  ex¬ 
presses  the  ratio  of  the  length  0V  to  the  intercept  0W  of  the  vector  v  =  0V  by  the 
covector  w  ,  w'here  W  means  the  intersection  of  0V  or  its  extension  with  the  termi¬ 
nal  hyperplane  of  the  covector  w  [S-RC8-9]  : 

v  w  =  v^w,  +  .  .  .  +  v\, ,  =  v°w  =  0V-.0W  .  a-10 

—  1  N  n 

When  the  vector  v  is  in  parallel  with  the  hypcrplane,  namely,  both  0  and  V  arc  in 
the  initial  hyperplane,  this  product  becomes  zero,  wrhich  is  the  case  of  t  =  0  in 
equation  a-4/2.  When  the  intersecting  point  W  is  on  the  extension  of  VO  =  -  v 
beyond  the  origin  0,  the  product  becomes  negative. 

Hence,  the  equation 

n  i  N 

w  x  =  w  x  +  .  . .  +  w  x  =  1  a-2 

n  1  N 
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expresses  a  hyperplane  passing  the  point  l/wn  of  the  n-th  axis,  and  xR  satisfying 

w  x11  =  t  ;  t  =  0,  1,  .  .  .  ,  T  a-4 

n  ... 

expresses  a  hyperplane  passing  through  the  point  t/w  of  n-th  axis,  in  parallel  to 
that  of  a- 2.  If  all  of  were  unity: 


w 

n 


=  1, 


n  =  1,  2, 


N 


a- 7 


the  resulting  equation 


12  N 

X  +X  +...+X  =1 


a- 8 


expresses  a  diagonal  hyperplane  passing  through  all  base  points  A  ,  A  ,  .  . .  ,  A 

1  u  IN 

and 


12  N 

X  +  X  +...+X  =t 


t  =  0,  1,  ...  ,  T 


a- 9 


expresses  diagonal  hyperplanes  intersecting  with  n-axis  at  the  length  t  hence  in 
parallel  to  that  of  a-8. 

If  the  base  vectors  are  expressed  as 


"a  =  0A  n  =  1,  2,  .  .  .  ,  N  , 
n  n 

base  covectors  as 


a-1 


a- 6 


then  there  holds 


— ■  m  -  m  ,  „  XT 

a  a  =1  ,  m,  n  =  1,  2,  .  .  .  ,  N  a-11 

n  — *  n 

where  1  m  means  the  unit  matrix  if  the  left  hand  side  a  a”1  is  regarded  as  a 
n  m  n 

matrix.  As  an  equation  in  the  index  notation,  1  is  identical  with  the  Kronecker 

symbol.  Unless  an  orthonormal  frame  is  defined  or  a  covector  w  is  substituted 
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by  a  vector  w  under  a  certain  geometric  correspondence,  an  "inner  product 
between  two  vectors"  and  hence  also  an  "angle"  arc  not  defined  and  the  Pythagorean 
theorem  is  not  applicable.  If  only  the  affine  frames  are  used,  the  space  is  called 
the  "affine  space." 

Since  an  "N-dimensional  space"  for  any  natural  number  N  is  none  other  than 
a  set  of  points  and  substantially  means  that  any  combination  of  the  numerical  value 
of  N  variables  is  mapped  to  a  single  point  in  this  set,  any  N-spacc  and  N-geometric 
object  can  be  drawn  on  paper  as  the  projected  figure  of  the  N-dimensional  figure  into 
the  two-dimensional  space  in  exactly  the  same  way  as  practiced  in  the  graphics  and 
descriptive  geometry.  Especially,  if  the  geometric  figures  are  discrete,  this  pro¬ 
jected  realization  in  2-space  can  be  proper  tw'o-dimensional  figures. 

N-space  (N  >  3)  is  already  often  used  among  engineers  [A-gemsf,  H-gs  17  2, 
183-98]  .  Physical  realization  of  N-space  by  means  of  electric  networks  w'as  ex¬ 
plained  in  1940  by  L.  Gutenmacher  in  G-kmuk  [see  also  K-NG,  S-MG-1,  -  2,  S-NG, 
C-RP  141,  etc.]  . 
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SECTION  III 


MULTIPLICATION  OF  THREE  NUMBERS 

The  geometrical  arrangement  of  this  case  is  exemplified  in  Figure  9,  but  the 
carrying  circuits  in  Figure  9  are  for  the  binary  case  which  is  explained  later.  In 
case  of 


n  =  A,  B, 

c 

1/3 

P  =  ABC 

2/3 

A  =  AjlO1  +  .  . 

.  +  A.101  +  .  .  . 

1 

+  A0.  V0.  A.  <10 

5.  1 

B  =  B  10J  +  .  . 
J 

.  +  B.10J  +  .  .  . 
3 

+  W°- V10 

5.  2 

c  =  c  ioK  + . 

XX 

k 

.  .  +  C.  10  +  .  ( 
k 

•+W°-Ck<10- 

5.  3 

The  product  P  is  obtained  by  the  carrying-over  process  from  the  quasi-decimal 
expression  Q  defined  by  equations  13  and  15  and  from  the  following  equation  14/3 
and  the  sum 


Q  =  2.  R 

t  i,  j  ;k  ijk 


t  =  0,  1,  ...  ,  T 


T  =  I  +  J  +  K 


10/3 


13 


14/3 


A.B.C. 
i  J  k 


over  all  the  values  of  i,  j,  and  k  satisfying 

i  +  j  +  k  =  t 


i  =  0,  1,  ...  ,  I 

j  =  0,  1,  . .  .  ,  J 
k  =  0,  1 . K  . 
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11/3 


12/3 
6.  1 
6.  2 


6.3 


I 

o 


o 

Ui 


I 


D 

a 

z 


X 
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Figure  9.  Multiplier  of  Three  Numbers  -  A,  B,  and  C 


If  any  three-dimensional  affine  frame  is  constructed  as  shown  in  the  three- 


dimensional  Figure  9  and  we 


1.  write  down  the  value  A.  at  all  points  of  the  coordinates: 
(x2,  x3)  =  (j  ,  k),  (j  =  0,  1,  .  .  .  ,  J  ;  k  =  0,  1,  .  .  .  ,  K) 


on  the  plane: 


x  =  i,  (i  =  0,  1,  . . .  ,  I) 


2.  write  down  the  value  Eh  at  all  points  of  the  coordinates: 
(x1,  x3)  =  (i,  k)  (i  =  0,  1 . I  ;  k  =  0,  1,  .  .  .  ,  K) 


on  the  plane: 


x  =  j,  (j  =  0,  1 . J) 


3.  write  down  the  value  C,  at  all  points  of  the  coordinates: 

rC 

(X1,  X2)  =  (i,  j)  (i  =  0,  1,  . . .  ,  I  ;  j  =  0,  1,  . . .  ,  J) 


on  the  plane: 


then 


x*  =  k,  (k  =  0,  1,  .  .  .  ,  K) 


4.  the  summands  of  lie  on  the  diagonal  plane 


a- 9/3: 


12  3 

X  +  X  +  X 

=  t 

t  =  0,  1,  .  .  . 

,  T 

T  =  I  +  J  +  K 


18.  1/3 


19.  1/3 


18.  2/3 


19.  2/3 


18.  3/3 


19.  3/3 


20/3 

20/3 
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From  these,  the  procedures  for  multiplication  are  established  as  follows: 


1. 


Place  A.  sheets  of  planes  at 

1 


1 

x  =  l,  i  =  0,  1,  .  . .  ,  I 


21. 1/3 


in  parallel  to  the  plane 


x  =  0 


namely,  the  initial  plane  of  a  covector  of  components 


(w1,  w2,  wg)  =  (1,  0,  0)  ; 


2.  Place  B.  sheets  at 
J 


x  =  j,  j  =  0,  1,  . . .  ,  J 


in  parallel  to  the  plane 


x  =  0  ; 


3.  Place  C,  sheets  at 
k 


x  =  k,  k  =  0,  1,  .  . .  ,  K 


22.  1/3 


23.  1/3 


21.  2/3 


22.  2/3 


21.  3/3 


in  parallel  to  the  plane 

3 

x  =  0  .  22.  3/3 

4.  Add  the  number  of  intersection  points  of  these  sets  of  plane  on  the 

diagonal  planes  of  equation  19/3  with  carry  from  over  to  Q  for  t  <  T  and  from 

O  to  P  or  further  to  P  „  etc.  Then  the  product  P  is  obtained.  The  total 
TT+1  T+2 

number  of  the  product  terms  A.B^C^  is  given  by 

M  =  (I  +  1)  (J  +  1)  (K  +  1)  .  7/3 

O 
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It  all  three  numbers  are  of  D  digits: 


1  +  1  =  J+1  =  K  +  1  =  D 


8/3 


then  there  holds 


M3 


9/3 
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SECTION  IV 


MULTIPLICATION  OF  N  NUMBERS 


In  case 


n  =  A,  B, . . .  ,  N  ; 


P  =  AB  . . .  n  . .  .  N  , 


In  n 

n  =  n  10  +...  +n  10  +...  +  n  ,  n  /  0,  n  <10 

Tn  .n  0  n  ’  .n 

I  l  I  l 

i11  =  0,  1,  ...  ,  In  ;  n  =  A,  B . N 


the  quasi-decimal  expression  Q  of  equation  15  is  defined  by  the  sum  : 


O  =  2  R 

t  .A  .n  .N  .A  .n  " '  .N 

i  i  1111 


10 


R  . .  =  A  .  .  .  n  .  .  .  N 

.A  .n  .N  .A  .n  .N 

1111  i  i 


11 


over  the  values  satisfying 


.A  .n  .N 

l  +...+1  +...+1  =  t  , 


12 


t  =  0,  1,  ...  ,  T  ; 


13 


m  tA  n  N 

T  =  I  +...+I  +...+I 


14 


1.  If  any  N-dimensional  affine  frame: 


V  a2 . a„ . “n 


a-1 


is  constructed,  either  electrically  or  in  any  other  manner  and 


22 


2.  the  value 


n  n  (n  A,  B,  .  .  .  ,  N  ;  in  =  0,  1 . 1°) 

i 

is  entered  at  all  points 


1  n  - 1  n  +  1  N  A 

(X  ,  .  .  .  ,  X  ,  X  ,  .  .  .  ,  X  )  =  (1  , 


.n  -  1  .n  +  1  ,N . 

1  ,i 
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on  the  hyperplane 


n 

x  = 


19 


3.  then  by  adding  the  products  at  first  on  the  diagonal  hyperplane  de¬ 
fined  by 

a’4-  a-7:  V“  -  t,  wn  -  1  (n  >  1.  2 . N) 

namely, 

1  2  n  N 

a-9:  x+x+...+x+...+x=t  20 


t  =  0.  1. 


20 


14: 


-i* 


+  In  + 


+  I 


N 
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Q  is  obtained  and  by  adding  these  Q  with  carry  over,  the  product  P  is  obtained, 
t  t 

Hence  the  rules  of  multiplication  are: 

1.  Construct  N-dimensional  affine  frame,  where  N  is  equal  to 
the  number  of  the  factors  A,  B,  .  .  .  ,  N  ; 

2.  Realize  n  (i11  =  0,  1,  . .  ,  ,  1°  ;  n  =  A,  B,  .  . ,  ,  N)  sheets  of 
hyperplanes  in  parallel  to  the  base  covector  a°  (n  =  1,  2,  .  . .  ,  N) 
at 


n 

x 


n  ,n  „  „  Tn, 

i  (l  =  0,  )  ; 


19.  6 
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3.  is  obtained  by  adding  the  number  of  intersections  of  these 
sets  of  hyperplanes  on  the  diagonal  hyperplanes  of  equation  19; 

4.  The  product  P  is  obtained  by  adding  these  Q^IO*  with  carry  over. 
If  the  requested  maximum  number  of  digits  is  D  for  all: 

1°  +  1  =  i  +  i  =  d 


T  =  NI  =  N(D  -  1) 
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by  installing 


mn  -  D 


N 


elements  of  multiplying  N  numbers,  all  M  multiplication  of 


R. .  =  A  .  .  .  n  . . .  N  ;  l  =0,1 . I 

.A  .n  .N  .A  .n  .N 

1111  l  l 


from 


R 


00 


=  A, 


N 


to 


R 


II  .  .  .  I 


.  N 


I 


can  be  simultaneously  done  spending  the  same  period  to  the  above  single 
multiplication. 

For  the  addition  of  these  products  to  Q  of  each  t,  instead  of  repeating 
the  addition  of  two  numbers  if  the  combinational  circuits  of  the  so-called  sym¬ 
metric  Boolean  functions  [C-SCLD236,  HiG-LODEC  83,  N-rc,  Sh-sarsc]  of 
Figure  10  are  used,  Q  can  be  obtained  in  a  single  step  independent  of  the  number 
of  its  digits.  Different  t  can  be  added  simultaneously. 
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Summands : 


Lo  =  s 1 ,  3 , 5 , 7 .  .  .  15 
1 1  =  S2, 3, 6, 7, 10, 11, 14, 15 
L2  =  S4,5,6, 7,12,13,14,15 
L3  =  S8,9,10 . 15 


O 

Sum:  five  digit  ^ 

binary  numbers  4 

Figure  10  Combinatorial  Adder 

When  the  number  of  the  summands  are  many,  e.  g.  ,  1024,  by  splitting  these 
into  8  groups  and  then  further  splitting  each  group  again  into  8  subgroups,  the 
addition  is  reduced  to  that  of  16  numbers  and  two  further  addition  of  eight  numbers. 
By  such  multiple-step  splitting,  the  combinational  circuits  of  the  type  of  Figure  10 
generally  can  be  kept  to  a  reasonable  size. 

Parallel  operation  of  the  addition  in  these  8  x  8  =  64  subgroups  reduces  the 
computing  time.  Series  operation  saves  the  adding  components  at  some  increase 
in  computing  time.  Because  the  purpose  of  this  paper  is  to  describe  the  basic 
ideas,  only  a  few  remarks  will  be  given  on  the  practical  realization. 

At  first,  the  error  correcting  and  detecting  facilities  can  be  introduced  by 
redundancy.  In  Figures  3  and  4,  if  nine  conductors  are  installed  instead  of  the 
bars,  and  magnetic  cores  are  linked  at  all  cross  points,  this  can  form  a  core 
matrix  and  the  multiplication  can  be  done  in  a  single  time  unit  by  impressing  the 


6  6  6  6 


binary  number 
input 

< 


lRn>0- 


25 


pulses  simultaneously  on  A.  and  B. .  The  conductors  for  each  decimal  digit  can 
be  reduced  to  five  as  explained  after  2:  equation  9/2.  The  decimal  value  of  each 
product  A.B^  can  be  obtained  (without  delay  caused  by  the  carrying  process)  by  thi 
use  of  the  combinatorial  circuits  of  the  type  shown  in  Figure  10  where  each  core 
of  9  x  9  =  81  or  5  x  5  =  25  forms  the  input.  The  addition  of  these  to  can  be  don 
again  by  the  combinatorial  circuits  of  type  Figure  10  without  carrying  delay. 


This  can  be  generalized  to  the  case  N,  where  the  magnetic  cores  link  to  the 
conductors  of  N-dimensional  lattice  circuit.  The  circuit  of  Figures  11,  12,  and 
13  is  a  combinatorial  circuit  whose  output  lines  P.  will  assume  directly  the  values 
that  represent  the  product  of  the  input  lines  A  and  B  in  the  binary  expression,  after 
a  suitable  logical  delay  for  propagation  of  voltage  levels  through  the  matrix.  For 
example: 


P  =  A  x  B  =  10110  x  01101  =  0100011110  ,  25/2 

The  highest  power  T  of  the  quasi-binary  product  Q  is 

T  =  I+  J=  4+  4  =  8;  26/2 

hence  by  possible  carry,  the  highest  power  S  of  the  product  P  is 

S  =  T  +  1  -  9  ,  27/2 


which  decimally  means 

P  =  A  x  B  =  22  x  18  =  286  28/2 

T  =  I+  J  =  1+1  =  2,  S  =  T  =  2  .  29/2 

In  Figures  12  and  13  the  boxes  in  general  position  have  three  inputs: 

i)  the  product,  namely  the  output  of  AND  gate  of  A^  and  By 

R..  =  A.B.  30/2 

ij  1  J 
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1  1 


o 

# 


Figure  II  Multiplier  of  Two  Binary  Numbers 


G  =0:  G .  „  G_  G.rt  G1A  G 

io  40  30  20  10  oo 


«-H  *  0 


Figure  12  Circuit  Construction 
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Figure  13  General  Cell 

ii)  input  from  the  upper  box:  G_ 

iii)  input  from  the  right  box:  H_ 

and  two  outputs: 

i)  adding  S  gate  to  lower  box: 

1,  3 

G.  =  R..  +  G..  +  H..  (mod.  2)  31/2 

l  - 1  j  +1  ij  ij  ij 

ii)  carrying  S  gate  to  left  box: 

l,  3 

H.  ,  .  =  G..H..  +  G..R.. +H..R..  (mod.  2)  32/2 

i+l  J  ij  ij  ij  ij  ij  ij 
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The  first  row  of  the  first  partial  product  has  no  input  from  above: 


and  no  output  to  the  left: 


Gi0  = 


33/2 


Hi+1  0  *  0  34/2 

and  further  no  input  from  right;  however,  this  condition  logically  results  from 
equations  30/2,  31/2,  and  32/2: 


Hi0  Gi  -  1  0Hi-l  0  +  Gi-l  0Ri-l  0  Hi  -  1  0Ri  -  1  0 


=  0-  0+  0-  R.  +  0  •  R.  ,  =0 

1-10  1-10 


The  first  right  skew  column  including  the  factor  has  no  input  from  right 


35/2 


H0J  -  0  • 


36/2 


The  input  of  the  first  left  skew  column  comes  from  the  adjacent  skew  upper 

right  horizontal  output: 


GT  .  =  Ht  „  . 

I J+l  1+1  J 


37/2 


For  j  =  0,  this  means 


GI  1  "  HI  +  1  0  =  ° 


38/2 


by  equation  31/2. 

The  first  J  digits  P  ,  P  . . .  ,  P  from  the  right  of  the  product  P  is  given  by 

0  1  J 

P.  =  G  .  .  .  ;  39/2 

J  "I  J+l 

further  digits  are  given  by 


i+J+1 


GiJ+l 


40/2 


29 


If  Y  is  the  time  required  to  form  the  logical  product  R_  from  A.  and 

and  Z  is  the  time  required  to  form  G.  „  .  ,  and  H.  ,  .  from  the  inputs 

l-l, j+1  i+l,3 

R„,  G„  and  H„  through  a  single  cell  of  the  matrix,  the  longest  delay  occurs  in 
the  series 


Hll’  «21'  H31-  H41’  °32'  H42-  H42'  G33'  «43'  °34'  H44'  H54  '  41/2 


hence  the  maximal  delay  U  for  the  entire  matrix  is  given  by 


U  =  Y  +  Z  |  (D  -  1)  +  2(D  -  2)  +  1 1  =  Y  +  Z(3D  -  4)  . 


42/2 


The  case  N  =  3  can  be  constructed  by  a  similar  circuitry,  but  with  some  additional 
gates  shown  by  the  small  circles  at  the  left  end  of  each  horizontal  plane  in  Fig¬ 
ure  9.  The  factors  of  the  same  radix  power  appear  in  the  same  vertical  diagonal 
planes  in  Figure  9: 


Rooo  VoVRioo  AiBoco’Roio  AoBico’ 


R„„,  =  AaBaC,,  ...  ,  R,,,  =  ARC  . 
001  001  444  444 


Because 


I  =  J  =  K  =  4 

5 

in  this  example,  all  of  the  A,  B,  and  C  are  smaller  than  2  ;  hence  the  product  P  is 
5  3  15 

smaller  than  (2  )  =  2  .  Therefore,  the  maximal  value  of  S  in  P  is  14: 

S  <  3(1  +  1)  -  1  =  14  . 
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SECTION  V 


CONCLUSION 


The  above  process  can  be  summarized  as  follows: 

1.  Simultaneous  expression  of  each  numeral  in  all  places  where- 
ever  it  is  used, 

2.  Mutual  separation  of  the  products  of  every  pair  or  set  of 
digits  by  removing  the  carryover  process  into  the  final 
adding  step, 

3.  Utilization  of  the  simple  "diagonal"  location  of  all  the  products 
of  pairs  or  sets  of  digits  of  the  same  radix  power,  and 

4.  Simultaneous  partial  addition  of  the  products  to  the  quasi-term 


by  the  combinational  circuits  instead  of  the  series  operation 
of  repeating  the  shift-and-add  with  carry  operations  of  two 
partial  products. 


By  assembling  these  geometrized  arithmetic  devices,  any  polynomial  or 
Taylor  series  of  finite  terms  can  be  realized  in  hardware,  which  will  facilitate 
and  expedite  the  numerical  calculation  of  roots  of  algebraic  equations. 


Satio  Okada 
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